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Abstract 

We study the unitary matrix model with a topological term. We call the 
topological term the theta term. In the symmetric model there is the phase 
transition between the strong and weak coupling regime at Ac = 2. If the Wilson 
term is bigger than the theta term, there is the strong-weak coupling phase 
transition at the same Ac. On the other hand, if the theta term is bigger than 
the Wilson term, there is only the strong coupling regime. So the topological 
phase transition disappears in this case. 
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1 Introduction 



Models of the symmetric unitary matrix are solved exactly in the double scaling 
limit, using orthogonal polynomials on a circle. The partition function is the 
form / dC/exp{-^tr([/ + ?7t)}, where U is an N X N unitary matrix. We call 
the model symmetric model. 

This unitary models has been studied in connection with the large- approx- 
imation to QCD in two dimensions. {one-plaquette model)^ Gross and Witten 
discovered the third-order phase transition between the weak and strong cou- 
pling regime at = 2. 

We consider the model which has the symmetric and anti-symmetric part. 
The symmetric part is the usual Wilson action. The anti-symmetric part be- 
comes the topological term. We can see the topological meaning of the theta 
term in the continuous limit, it gives rise to a phase transition at 6 = n, if the 
Wilson term is bigger than the theta term. Q We call this phase transition the 
topological phase transition and the model which has only the theta term the 
anti-symmetric model. 

From the view point of the integrable system this model can be embed- 
ded in the two-dimensional Toda hierarchy with the conjugate structure. [^,0 
This Toda equation is split into the modified Volterra (MV) equation and the 
discrete nonlinear Schrodinger (DNLS) equation. MV and DNLS equations cor- 
respond to the symmetric and anti-symmetric model respectively. Coupling the 
Toda equation and the string equation we can derive the third Painlevc (P III) 
equation. We use P III to study the phase structure. 

This letter is organized as follows. In the section 2 we introduce the unitary 
matrix model with a topological term. We consider the two cases: 

(i) the Wilson term is bigger than the theta term, 

(ii) the theta term is bigger than the theta term. 

In the section 3 we consider the phase transition in the case (i) using the 
P III. In the section 4 we study the phase structure in the case (ii). The last 
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section is devoted to the concluding remarks. 



2 Wilson term and theta term 

We consider the partition function of the unitary matrix model with a topolog- 
ical term. We consider the unitary matrix model 

f N 

Zn= dUeM-jV{U)), (2.1) 

where V{U) is a potential 

V{U) =tiU + t-iU''. (2.2) 

U is the gauge group U{N). 

We divide the potential into the symmetric and the anti-symmetric part, 

V{U) = t+s^ + t_S0, (2.3) 

where 

ti + t-i ti — t-i 

*+ = ^— ' *- = ^— (2.4) 

Sw is the symmetric part, the usual Wilson action 

= l(trC/ + trZ/t). (2.5a) 



Here we choose 



= i(trC/-trC/t), (2.5b) 



for the theta term, the anti-symmetric part. 

Hereafter we call next reduced models the symmetric model and the anti- 
symmetric model: 

t\ = t-i — t_ = 0, (symmetric model) 

ti = —t-i=t-, = 0. (anti — symmetric model) 
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Here we consider the two cases 

(i)i<+i > m-\>\t+\- (2.6) 

We parameterize ti and t-i by e: 

(i) ti = -e^ t-i = -e~% 

(ii) ^l = -e^ i_i = e-^ (2.7) 
The measure dU may be written as 

dU^l[^A{a)A{a). (2.8) 

m 

Here the eigenvalues of U are {exp(iai), exp(ia2), • ■ ■ , exp(iQ;jv)} and AA is 
the Jacobian for the change of variables, 

A(a) = detjfee'"^(^-'=\ 

A(a) = detjfce-'"^-^^-'^'). (2.9) 

Then we obtain the partition function in the case (i): 

Zn = const.detjfee'("^'+''')/_j+fe(iV/A) = const.detjfc/_j+fc(A^/A). (2.10) 

Here I-j+k is the modified Bessel function of order — j + k. In the same way we 
can calculate the partition hmction in the case (ii): 



const. detjfce'^ J_j+fe(iV/ A) const.detjfe J_j+fc(^/A). (2.11) 



Here J-j+k is the Bessel function of order —j + k. Notice that (2.1C) and (2.11) 
do not depend on e. 

It is well known that the partition function Zn of the unitary matrix model 
can be presented as a product of norms of the biorthogonal polynomial system. 
Namely, let us introduce a scalar product of the form 

<A,B>=i ^ eM-V{z)}A{z)B{z-'), (2.12) 



where 

V{z) ^tiz + t-iz^\ (2.13) 

Let us define the system of tlic polynomials biorthogaonal with respect to this 
scalar product 

<^n,^l>=KSnk- (2.14) 

Then, the partition function is equal to the product of /i„'s: 

Zn=1[ hk, Zo = 1. (2.15) 

fc=0 

The polynomials are normalized as follows (Note that the asterisk '*' does not 
mean the complex conjugation): 

$„ = 2" + --- + 5„-i, + + 5-1-5*1 = 1. (2.16) 

Now it is easy to show that these polynomials satisfy the following recurrent 
relations, 

$„+i(z) = z$„(z) + 5„z"$,*(z~i), 

= z-'Ki^-') + S*„z.Mz), (2.17) 



and 

hn+l 



From ( ^.14 ) we can obtain the string equations 



= l-S^Sl. (2.18) 



{n + l)Sn = {t-iSn+i + tiSn-i){l ~ SnSl), (2.19a) 

(n + l)Sl = ihS*+, + t-i5*_i)(l - 5„5*). (2.19b) 

In the unitary matrix model there is a conjugate relation: 

tlSnS^_1 — t-lS^Sn-l- (2.20) 
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Here we define a„ 



an = l- -SnS: = (2.21) 

/1 77, 



From ( p.20D a„ are functions of the radial coordinate 

X = iit_i, (2.22) 

only. 

a„ satisfies the next Painleve V with Sy =0:0 

9a;^ 2 a,i — 1 a„ 9a:: a; 9a:: 

2 ^ , (n+l)"a„-l 
a„(a„-l)H — ^ . (2.23) 

3 Phase structure in the case (i) 



The partition function dose not depend on e from (2.10). This can be seen from 



the radial coordinate ( ^.22 ). From these results in the large- iV limit the phase 



structure in the case (i) is the same as the symmetric model. To study the 
strong- weak coupling phase transition we use ( 2.23| ) in a; — > oo. We rewrite 
(|^) a second order ODE for S„ 



where i+ = N/X. This equation can be obtained directly from coupling the 
string equation and the modified Volterra equation. Q 

In particular when we consider the strong coupling regime we only need to 
solve 

dSn , 1 dSn An + lf 



dt\ t+ dt+ ^ t] 



4]5„ = 0. (3.2) 



This is the Bessel equation. Then setting n = N, we can obtain in the large- iV 
limit 

Sn = JN{2t+) + Oij^^) ''-^ JNi2N/X), (3.3) 
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where Jjv is the standard Bessel function. As a consequence, we find 



(3.4) 



A2 ° 2 
is especially appropriate for a discussion of the weak coupling regime. 
One may consider an 1/N expansion for Sat:!^ 

1 A3 



S 



N 



1 - - - — — (1 - + 0( — ) 

2 N^128^ 2^ ^^iV4^" 



(3.5) 



We notice the typical phase transition behavior of Sn- The critical point Ac = 2 
is independent of the parameter e. 

4 Phase structure in the case (ii) 

In the large iV-limit to study the phase structure we use the string equations 
( ^.19a ) and ( ^.19b| ). Setting n = N, there is a critical point at A^ = 2sinhe for 
the roots of the differential equation are degenerate. It seems that in the limit 
e — > (the Wilson term vanishes) the weak coupling regime disappears. But 
as seeing in the previous section this is not correct. From the radial coordinate 



or (2.11) the critical point docs not depend on e. The phase structure in the 
case (ii) is the same as the anti-symmetric model. To see the phase structure 
we study ( p3| ) in a; — CX3. Then the difference betweeen the case (i) and (ii) 



is the sign of the third term of RHS of (2.23). We rewrite (2.23) a second order 
ODE for Sn- 



1 - ST at- 



1 dSn , {n + If Sn 



1 s?, 



and 



1 dSn , [n + If S„ 



1 + SI 



+ 45„(1-S2), 
{n = odd) (4.1a) 

+ 4Sn{l + Sl), 
(n = evenjA.lh) 
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using 

Sn = S*, a„ = 1 ~ Sl, {n = odd), 

5„ = -S*„, a„ = l + 5'2, (n = even). (4.2) 

Notice t- = 2N/X. This equation can be obtained directly from the couphng 
the string equations and the discrete nonhnear Schrodinger equation, (t) When 
we consider the strong couphng regime we only need to solve 
9Sn 1 dS^ (n + ir _ 

This is the modified Bessel equation. Then setting n — N, we can obtain in the 
large- limit 

Sn = KN{2t^) + Oij^^) ""-^ Kn{2N/X), (4.4) 

where Kn is the second kind modified Bessel function. As a consequence, we 
obtain 

4 
A2 



Sn -e^pN\jl + -^. (4.5) 



On the other hand we can not do an 1/A^ expansion in (4.1a) and (4. It), since 
there is not the weak coupling regime. Then there is not the strong-weak cou- 
pling phase transition in the case (ii). 

In the previous letter, Q adding the term I log U we have shown that there 
is the phase transition at = tt in the case (i). We call the phase transition the 
topological phase transition. This phase transition can be seen in the weak cou- 
pling regime. Q So in the case (ii) there is not the topological phase transition, 
too. 

5 Concluding remarks 

We study the unitary matrix model with a topological term. This model con- 
tains the Wilson term and the theta term. The Wilson term is the symmetric 
and the theta is the anti-symmetric part. We call the model which has only the 
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Wilson term the symmetric model and the model which has only the theta term 
the anti-symmetric model. It is well known that in the symmetric model there 
is the strong- weak coupling phase transition at Ac = 2. If the Wilson term is 
bigger than the theta term, there is the strong-weak coupling phase transition at 
the same Ac- Adding the term llogU, there is the topological phase transition 
ed, 6 = TT. On the other hand if the theta term is bigger than the Wilson term, 
there is no strong-weak coupling phase transition. Since there is only the strong 
coupling regime, the topological phase transition also disappears. 
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